Analysis Of Bacterial Genomes Research Articles

Modeling the mosaic structure of bacterial genomes to infer their evolutionary history

The chronology and phylogeny of bacterial evolution are difficult to reconstruct due to a scarce fossil record. The analysis of bacterial genomes remains challenging because of large sequence divergence, the plasticity of bacterial genomes due to frequent gene loss, horizontal gene transfer, and differences in selective pressure from one locus to another. Therefore, taking advantage of the rich and rapidly accumulating genomic data requires accurate modeling of genome evolution. An important technical consideration is that loci with high effective mutation rates may diverge beyond the detection limit of the alignment algorithms used, biasing the genome-wide divergence estimates toward smaller divergences. In this article, we propose a novel method to gain insight into bacterial evolution based on statistical properties of genome comparisons. We find that the length distribution of sequence matches is shaped by the effective mutation rates of different loci, by the horizontal transfers, and by the aligner sensitivity. Based on these inputs, we build a model and show that it accounts for the empirically observed distributions, taking the Enterobacteriaceae family as an example. Our method allows to distinguish segments of vertical and horizontal origins and to estimate the time divergence and exchange rate between any pair of taxa from genome-wide alignments. Based on the estimated time divergences, we construct a time-calibrated phylogenetic tree to demonstrate the accuracy of the method.

Bacterial isolation and genome analysis of a novel Klebsiella quasipneumoniae phage in southwest China’s karst area

BackgroundSouthwest China is one of the largest karst regions in the world. Karst environment is relatively fragile and vulnerable to human activities. Due to the discharge of sewage and domestic garbage, the karst system may be polluted by pathogenic bacteria. The detection of bacterial distribution and identification of phage capable of infecting them is an important approach for environmental assessment and resource acquisition.MethodsBacteria and phages were isolated from karst water in southwest China using the plate scribing and double plate method, respectively. Isolated phage was defined by transmission electron microscopy, one-step growth curve and optimal multiplicity of infection (MOI). Genomic sequencing, phylogenetic analysis, comparative genomic and proteomic analysis were performed.ResultsA Klebsiella quasipneumoniae phage was isolated from 32 isolates and named KL01. KL01 is morphologically identified as Caudoviricetes with an optimal MOI of 0.1, an incubation period of 10 min, and a lysis period of 60 min. The genome length of KL01 is about 45 kb, the GC content is 42.5%, and it contains 59 open reading frames. The highest average nucleotide similarity between KL01 and a known Klebsiella phage 6939 was 83.04%.ConclusionsKL01 is a novel phage, belonging to the Autophagoviridae, which has strong lytic ability. This study indicates that there were not only some potential potentially pathogenic bacteria in the karst environment, but also phage resources for exploration and application.

Functional analyses of bacterial genomes found in Symbiodiniaceae genome assemblies.

Bacterial-algal interactions strongly influence marine ecosystems. Bacterial communities in cultured dinoflagellates of the family Symbiodiniaceae have been characterized by metagenomics. However, little is known about whole-genome analysis of marine bacteria associated with these dinoflagellates. We performed in silico analysis of four bacterial genomes from cultures of four dinoflagellates of the genera Symbiodinium, Breviolum, Cladocopium and Durusdinium. Comparative analysis showed that the former three contain the alphaproteobacterial family Parvibaculaceae and that the Durusdinium culture includes the family Sphingomonadaceae. There were no large genomic reductions in the alphaproteobacteria with genome sizes of 2.9-3.9 Mb, implying they are not obligate intracellular bacteria. Genomic annotations of three Parvibaculaceae detected the gene for diacetylchitobiose deacetylase (Dac), which may be involved in the degradation of dinoflagellate cell surfaces. They also had metabolic genes for dissimilatory nitrate reduction to ammonium (DNRA) in the nitrogen (N) cycle and cobalamin (vitamin B12 ) biosynthetic genes in the salvage pathway. Those three characters were not found in the Sphingomonadaceae genome. Predicted biosynthetic gene clusters for secondary metabolites indicated that the Parvibaculaceae likely produce the same secondary metabolites. Our study suggests that the Parvibaculaceae is a major resident of Symbiodiniaceae cultures with antibiotics.

Bacterial community and genome analysis of cytoplasmic incompatibility-inducing Wolbachia in American serpentine leafminer, Liriomyza trifolii.

Liriomyza trifolii, an agricultural pest, is occasionally infected by Wolbachia. A Wolbachia strain present in Liriomyza trifolii is associated with cytoplasmic incompatibility (CI) effects, leading to the death of embryos resulting from incompatible crosses between antibiotic-treated or naturally Wolbachia-free strain females and Wolbachia-infected males. In this study, high-throughput sequencing of hypervariable rRNA genes was employed to characterize the bacterial community in Wolbachia-infected L. trifolii without antibiotic treatment. The analysis revealed that Wolbachia dominates the bacterial community in L. trifolii, with minor presence of Acinetobacter, Pseudomonas, and Limnobacter. To elucidate the genetic basis of the CI phenotype, metagenomic sequencing was also conducted to assemble the genome of the Wolbachia strain. The draft-genome of the Wolbachia strain wLtri was 1.35 Mbp with 34% GC content and contained 1,487 predicted genes. Notably, within the wLtri genome, there are three distinct types of cytoplasmic incompatibility factor (cif) genes: Type I, Type III, and Type V cifA;B. These genes are likely responsible for inducing the strong cytoplasmic incompatibility observed in L. trifolii.

Mosquitoes provide a transmission route between possums and humans for Buruli ulcer in southeastern Australia

Buruli ulcer, a chronic subcutaneous infection caused by Mycobacterium ulcerans, is increasing in prevalence in southeastern Australia. Possums are a local wildlife reservoir for M. ulcerans and, although mosquitoes have been implicated in transmission, it remains unclear how humans acquire infection. We conducted extensive field survey analyses of M. ulcerans prevalence among mosquitoes in the Mornington Peninsula region of southeastern Australia. PCR screening of trapped mosquitoes revealed a significant association between M. ulcerans and Aedes notoscriptus. Spatial scanning statistics revealed overlap between clusters of M. ulcerans-positive Ae. notoscriptus, M. ulcerans-positive possum excreta and Buruli ulcer cases, and metabarcoding analyses showed individual mosquitoes had fed on humans and possums. Bacterial genomic analysis confirmed shared single-nucleotide-polymorphism profiles for M. ulcerans detected in mosquitoes, possum excreta and humans. These findings indicate Ae. notoscriptus probably transmit M. ulcerans in southeastern Australia and highlight mosquito control as a Buruli ulcer prevention measure.

Animal Microbiomes as a Source of Novel Antibiotic-Producing Strains.

Natural compounds continue to serve as the most fruitful source of new antimicrobials. Analysis of bacterial genomes have revealed that the biosynthetic potential of antibiotic producers by far exceeds the number of already discovered structures. However, due to the repeated discovery of known substances, it has become necessary to change both approaches to the search for antibiotics and the sources of producer strains. The pressure of natural selection and the diversity of interactions in symbiotic communities make animal microbiomes promising sources of novel substances. Here, microorganisms associated with various animals were examined in terms of their antimicrobial agents. The application of alternative cultivation techniques, ultrahigh-throughput screening, and genomic analysis facilitated the investigation of compounds produced by unique representatives of the animal microbiota. We believe that new strategies of antipathogen defense will be discovered by precisely studying cell-cell and host-microbe interactions in microbiomes in the wild.

Genome-based reclassification of Kitasatospora niigatensis as a later heterotypic synonym of Kitasatospora cineracea Tajima et al. (2001).

The present study used genome-based approaches to investigate the taxonomic relationship between Kitasatospora cineracea DSM 44780T and Kitasatospora niigatensis DSM 44781T, two species that were previously described by Tajima et al. (Int J Syst Evol Microbiol 51:1765-1771, 2001). The digital DNA-DNA hybridization (dDDH), average amino acid identity (AAI), and average nucleotide identity (ANI) values between the genomes of the two type strains were 90.3, 98.7, and 99.1%, respectively. These values exceeded the established thresholds of 70% (dDDH) and 95-96% (ANI and AAI) for bacterial species delineation, suggesting that K. cineracea and K. niigatensis should share the same taxonomic position. Furthermore, our analysis using the 'Bacterial Pan Genome Analysis' (BPGA) pipeline and the Maximum Likelihood core-genes tree inferred using FastTree2 consistently demonstrated that K. cineracea DSM 44780T and K. niigatensis DSM 44781T are closely related, as indicated by the clustering of these strains in the core-genes phylogenomic tree. Based on these findings, we propose that K. niigatensis should be considered a later heterotypic synonym of K. cineracea.

Accurate and fast graph-based pangenome annotation and clustering with ggCaller.

Bacterial genomes differ in both gene content and sequence mutations, which underlie extensive phenotypic diversity, including variation in susceptibility to antimicrobials or vaccine-induced immunity. To identify and quantify important variants, all genes within a population must be predicted, functionally annotated, and clustered, representing the "pangenome." Despite the volume of genome data available, gene prediction and annotation are currently conducted in isolation on individual genomes, which is computationally inefficient and frequently inconsistent across genomes. Here, we introduce the open-source software graph-gene-caller (ggCaller). ggCaller combines gene prediction, functional annotation, and clustering into a single workflow using population-wide de Bruijn graphs, removing redundancy in gene annotation and resulting in more accurate gene predictions and orthologue clustering. We applied ggCaller to simulated and real-world bacterial data sets containing hundreds or thousands of genomes, comparing it to current state-of-the-art tools. ggCaller has considerable speed-ups with equivalent or greater accuracy, particularly with data sets containing complex sources of error, such as assembly contamination or fragmentation. ggCaller is also an important extension to bacterial genome-wide association studies, enabling querying of annotated graphs for functional analyses. We highlight this application by functionally annotating DNA sequences with significant associations to tetracycline and macrolide resistance in Streptococcus pneumoniae, identifying key resistance determinants that were missed when using only a single reference genome. ggCaller is a novel bacterial genome analysis tool with applications in bacterial evolution and epidemiology.

Nitrogen removal capability and mechanism of a novel heterotrophic nitrification–aerobic denitrification bacterium Halomonas sp. DN3

Recently, the functional microorganisms capable of eliminating nitrogenous waste have been applied in mariculture systems. As a potential candidate for treating mariculture wastewater, strain DN3 eliminated 100% of ammonia and nitrate and 96.61%–100% of nitrite within 72 h, when single nitrogen sources at concentrations of 0–50 mg/L. Strain DN3 also exhibited the efficient removal performance of mixed-form nitrogen (ammonia, nitrate, and nitrite) at salinity 30 ‰, C/N ratio 20, and 180 rpm. The nitrogen assimilation pathway dominated inorganic nitrogen metabolism, with less nitrogen (14.23%-25.02% of TN) lost into the air via nitrification and denitrification, based on nitrogen balance analysis. Moreover, the bacterial nitrification pathway was explored by enzymatic assays and inhibition assays. These complex nitrogen assimilation and dissimilation processes were further revealed by bacterial genome analysis. These results provide important insight into nitrogen metabolism of Halomonas sp. and theoretical support for treating mariculture wastewater with strain DN3.

BacSeq: A User-Friendly Automated Pipeline for Whole-Genome Sequence Analysis of Bacterial Genomes.

Whole-genome sequencing (WGS) of bacterial pathogens is widely conducted in microbiological, medical, and clinical research to explore genetic insights that could impact clinical treatment and molecular epidemiology. However, analyzing WGS data of bacteria can pose challenges for microbiologists, clinicians, and researchers, as it requires the application of several bioinformatics pipelines to extract genetic information from raw data. In this paper, we present BacSeq, an automated bioinformatic pipeline for the analysis of next-generation sequencing data of bacterial genomes. BacSeq enables the assembly, annotation, and identification of crucial genes responsible for multidrug resistance, virulence factors, and plasmids. Additionally, the pipeline integrates comparative analysis among isolates, offering phylogenetic tree analysis and identification of single-nucleotide polymorphisms (SNPs). To facilitate easy analysis in a single step and support the processing of multiple isolates, BacSeq provides a graphical user interface (GUI) based on the JAVA platform. It is designed to cater to users without extensive bioinformatics skills.

In Vitro Activity of Cefiderocol on Multiresistant Bacterial Strains and Genomic Analysis of Two Cefiderocol Resistant Strains.

Cefiderocol is a new siderophore cephalosporin that is effective against multidrug-resistant Gram-negative bacteria, including carbapenem-resistant strains. The aim of this study was to evaluate the activity of this new antimicrobial agent against a collection of pathogens using broth microdilution assays and to analyze the possible mechanism of cefiderocol resistance in two resistant Klebsiella pneumoniae isolates. One hundred and ten isolates were tested, comprising 67 Enterobacterales, two Acinetobacter baumannii, one Achromobacter xylosoxidans, 33 Pseudomonas aeruginosa and seven Stenotrophomonas maltophilia. Cefiderocol showed good in vitro activity, with an MIC < 2 μg/mL, and was able to inhibit 94% of the tested isolates. We observed a resistance rate of 6%. The resistant isolates consisted of six Klebsiella pneumoniae and one Escherichia coli, leading to a resistance rate of 10.4% among the Enterobacterales. Whole-genome sequencing analysis was performed on two cefiderocol-resistant Klebsiella pneumoniae isolates to investigate the possible mutations responsible for the observed resistance. Both strains belonged to ST383 and harbored different resistant and virulence genes. The analysis of genes involved in iron uptake and transport showed the presence of different mutations located in fhuA, fepA, iutA, cirA, sitC, apbC, fepG, fepC, fetB, yicI, yicJ, and yicL. Furthermore, for the first time, to the best of our knowledge, we described two Klebsiella pneumoniae isolates that synthesize a truncated fecA protein due to the transition from G to A, leading to a premature stop codon in the amino acid position 569, and a TonB protein carrying a 4-amino acid insertion (PKPK) after Lysine 103. In conclusion, our data show that cefiderocol is an effective drug against multidrug-resistant Gram-negative bacteria. However, the higher resistance rate observed in Enterobacterales underlines the need for active surveillance to limit the spread of these pathogens and to avoid the risks associated with the emergence of resistance to new drugs.

Inter- and intraspecific phytochemical variation correlate with epiphytic flower and leaf bacterial communities.

Microbes associated with flowers and leaves affect plant health and fitness and modify the chemical phenotypes of plants with consequences for interactions of plants with their environment. However, the drivers of bacterial communities colonizing above-ground parts of grassland plants in the field remain largely unknown. We therefore examined the relationships between phytochemistry and the epiphytic bacterial community composition of flowers and leaves ofRanunculus acris and Trifolium pratense. On 252 plant individuals, we characterized primary and specialized metabolites, that is, surface sugars, volatile organic compounds (VOCs), and metabolic fingerprints, as well as epiphytic flower and leaf bacterial communities. The genomic potential of bacterial colonizers concerning metabolic capacities was assessed using bacterial reference genomes. Phytochemical composition displayed pronounced variation within and between plant species and organs, which explained part of the variation in bacterial community composition. Correlation network analysis suggests strain-specific correlations with metabolites. Analysis of bacterial reference genomes revealed taxon-specific metabolic capabilities that corresponded with genes involved in glycolysis and adaptation to osmotic stress. Our results show relationships between phytochemistry and the flower and leaf bacterial microbiomes suggesting that plants provide chemical niches for distinct bacterial communities. In turn, bacteria may induce alterations in the plants' chemical phenotype. Thus, our study may stimulate further research on the mechanisms of trait-based community assembly in epiphytic bacteria.

Determination of prebiotic utilisation capability of potential probiotic Bacillus velezensis FS26 through in silico and in vitro approaches

Prebiotics are substrates that can increase bacterial growth while also benefiting the host. The presence of enzymes responsible for hydrolysing the glycosidic bond in the prebiotic structure is critical for evaluating the prebiotic's use by the bacterial strain. In silico analysis of bacterial genomes is a method for predicting the hydrolysing enzymes that cleave the glycosidic link in prebiotic structures. Furthermore, in vitro validation must be performed to further analyse and evaluate the bacterial ability to utilise prebiotics. Thus, this study aims to investigate the ability of potential probiotic Bacillus velezensis FS26 to utilise prebiotics using in silico and in vitro techniques. B. velezensis FS26 genome analysis has shown the presence of glycosidase hydrolase enzymatic family that are capable of hydrolysing α- and β-glycosidic bonds in polysaccharide and oligosaccharide structures. The enzymes responsible for prebiotic utilisation were found in the bacterial genomes including α-galactosidase, fructofuranosidase, levansucrase and β-galactosidase. Furthermore, an in vitro study employing the phenol red carbohydrates assay revealed that B. velezensis FS26 can utilise the prebiotics lactulose, raffinose, and inulin, with comparable results to the probiotic Lactiplantibacillus plantarum ATCC 8014. Furthermore, the addition of prebiotics raffinose and inulin in the nutrient broth significantly increased the growth of B. velezensis FS26 cultures in contrast with prebiotic lactulose and nutrient broth without prebiotics. The agreement of in silico and in vitro results indicates the promising capability of B. velezensis FS26 to utilise prebiotics in aquaculture, which is an important sector for maintaining global food security.

Ecological divergence of syntopic marine bacterial species is shaped by gene content and expression

Identifying mechanisms by which bacterial species evolve and maintain genomic diversity is particularly challenging for the uncultured lineages that dominate the surface ocean. A longitudinal analysis of bacterial genes, genomes, and transcripts during a coastal phytoplankton bloom revealed two co-occurring, highly related Rhodobacteraceae species from the deeply branching and uncultured NAC11-7 lineage. These have identical 16S rRNA gene amplicon sequences, yet their genome contents assembled from metagenomes and single cells indicate species-level divergence. Moreover, shifts in relative dominance of the species during dynamic bloom conditions over 7 weeks confirmed the syntopic species’ divergent responses to the same microenvironment at the same time. Genes unique to each species and genes shared but divergent in per-cell inventories of mRNAs accounted for 5% of the species’ pangenome content. These analyses uncover physiological and ecological features that differentiate the species, including capacities for organic carbon utilization, attributes of the cell surface, metal requirements, and vitamin biosynthesis. Such insights into the coexistence of highly related and ecologically similar bacterial species in their shared natural habitat are rare.

Analysis of multipartite bacterial genomes using alignment free and alignment-based pipelines

Since the discovery of second chromosome in Rhodobacter sphaeroides 2.4.1 in 1989, multipartite genomes have been reported in over three hundred bacterial species under nine different phyla. This has shattered the unipartite (single chromosome) genome dogma in bacteria. Since then, many questions on various aspects of multipartite genomes in bacteria have been addressed. However, our understanding of how multipartite genomes emerge and evolve is still lacking. Importantly, the knowledge of genetic factors underlying the differences in multipartite and single-chromosome genomes is lacking. In this work, we have performed comparative evolutionary and functional genomics analyses to identify molecular factors that discriminate multipartite from unipartite bacteria, with the goal to decipher taxon-specific factors, and those that are prevalent across the taxa, underlying these traits. We assessed the roles of evolutionary mechanisms, specifically gene gain, in driving the divergence of bacteria with single and multiple chromosomes. In addition, we performed functional genomic analysis to garner support for our findings from comparative evolutionary analysis. We found genes such as those encoding conserved hypothetical proteins in Deinococcus radiodurans R1, and putative phage phi-C31 gp36 major capsid like and hypothetical proteins in Rhodobacter sphaeroides 2.4.1, which are located on accessory chromosomes in these bacteria but were not found in the inferred ancestral sequences, and on the primary chromosomes, as well as were not found in their closest relatives with single chromosome within the same clade. Our study shines a new light on the potential roles of the secondary chromosomes in helping bacteria with multipartite genomes to adapt to specialized environments or growth conditions.

Dietary fiber chemical structure determined gut microbiota dynamics

AbstractPrecision modulation of gut microbiota requires elucidation of the relation between dietary fiber intake and gut microbe dynamics. However, current studies on this aspect are few due to many technical limitations. Here, we used Caenorhabditis elegans to minimize the complicated host–microbial factors and to find the relation between dietary fiber chemical structures and gut microbiota dynamics. The Allium schoenoprasum polysaccharide (AssP) structure was elucidated and used as the complex dietary fiber against the simple fiber inulin. In vitro bacterial growth and genome analysis indicated that AssP supports bacterial growth better than inulin, while in vivo gut microbiota analysis of C. elegans fed with AssP showed that microbiota richness increased significantly compared with those fed with inulin. It is concluded that the more complex the dietary fiber chemical structure, the more gut bacteria growth it supports. Together with the community bacterial interactions that alter their abundances in vivo, these factors regulate gut microbiota synergistically.

Integrated Pangenome Analysis and Pharmacophore Modeling Revealed Potential Novel Inhibitors against Enterobacter xiangfangensis.

Enterobacter xiangfangensis is a novel, multidrug-resistant pathogen belonging to the Enterobacter genus and has the ability to acquire resistance to multiple antibiotic classes. However, there is currently no registered E. xiangfangensis drug on the market that has been shown to be effective. Hence, there is an urgent need to identify novel therapeutic targets and effective treatments for E. xiangfangensis. In the current study, a bacterial pan genome analysis and subtractive proteomics approach was employed to the core proteomes of six strains of E. xiangfangensis using several bioinformatic tools, software, and servers. However, 2611 nonredundant proteins were predicted from the 21,720 core proteins of core proteome. Out of 2611 nonredundant proteins, 372 were obtained from Geptop2.0 as essential proteins. After the subtractive proteomics and subcellular localization analysis, only 133 proteins were found in cytoplasm. All cytoplasmic proteins were examined using BLASTp against the virulence factor database, which classifies 20 therapeutic targets as virulent. Out of these 20, 3 cytoplasmic proteins: ferric iron uptake transcriptional regulator (FUR), UDP-2,3diacylglucosamine diphosphatase (UDP), and lipid-A-disaccharide synthase (lpxB) were chosen as potential drug targets. These drug targets are important for bacterial survival, virulence, and growth and could be used as therapeutic targets. More than 2500 plant chemicals were used to molecularly dock these proteins. Furthermore, the lowest-binding energetic docked compounds were found. The top five hit compounds, Adenine, Mollugin, Xanthohumol C, Sakuranetin, and Toosendanin demonstrated optimum binding against all three target proteins. Furthermore, molecular dynamics simulations and MM/GBSA analyses validated the stability of ligand-protein complexes and revealed that these compounds could serve as potential E. xiangfangensis replication inhibitors. Consequently, this study marks a significant step forward in the creation of new and powerful drugs against E. xiangfangensis. Future studies should validate these targets experimentally to prove their function in E. xiangfangensis survival and virulence.

Single-cell genomics reveals population structures from in vitro evolutionary studies of Salmonella.

Single-cell DNA sequencing has the potential to reveal detailed hierarchical structures in evolving populations of cells. Single cell approaches are increasingly used to study clonal evolution in human ageing and cancer but have not yet been deployed to study evolving clonal microbial populations. Here, we present an approach for single bacterial genomic analysis for in vitro evolution experiments using FACS isolation of individual bacteria followed by whole-genome amplification and sequencing. We apply this to the experimental evolution of a hypermutator strain of Salmonella in response to antibiotic stress (ciprofloxacin). By analysing sequence polymorphisms in individual cells from populations we identified the presence and prevalence of sub-populations which have acquired polymorphisms in genes previously demonstrated to be associated with ciprofloxacin susceptibility. We were also able to identify that the population exposed to antibiotic stress was able to develop resistance whilst maintaining diversity. This population structure could not be resolved from bulk sequence data, and our results show how high-throughput single-cell sequencing can enhance experimental studies of bacterial evolution.

Small wards in the ICU: a favorable measure for controlling the transmission of carbapenem-resistant Klebsiella pneumoniae

PurposeCarbapenem-resistant Klebsiella pneumoniae (CRKP) is one of the leading causes of healthcare-associated infections (HAIs) and is particularly pervasive in intensive care units (ICUs). This study takes ICU layout as the research object, and integrates clinical data and bacterial genome analysis to clarify the role of separate, small wards within the ICU in controlling the transmission of CRKP.MethodsThis study prospectively observed the carriage and spread of CRKP from a long-term in-hospital patient (hereafter called the Patient) colonized with CRKP in the gut and located in a separate, small ward within the ICU. The study also retrospectively investigated CRKP-HAIs in the same ICU. The relationship and transmission between CRKP isolates from the Patient and HAI events in the ICU were explored with comparative genomics.ResultsIn this study, 65 CRKP-HAI cases occurred during the investigation period. Seven CRKP-HAI outbreaks were also observed. A total of 95 nonrepetitive CRKP isolates were collected, including 32 strains from the Patient in the separate small ward. Phylogenetic analysis based on core genome single-nucleotide polymorphism (cgSNP) showed that there were five possible CRKP clonal transmission events and two clonal outbreaks (A1, A2) during the study. CRKP strains from the Patient did not cause CRKP between-patient transmission or outbreaks in the ICU during the 5-year study period.ConclusionThe presence of a long-term hospitalized patient carrying CRKP and positioned in a separate, small ward did not lead to CRKP transmission or infection outbreaks in the ICU. Combining a small-ward ICU layout with normative HAI control measures for multidrug-resistant pathogen infection was effective in reducing CRKP transmission.

Kmer2vec: A Novel Method for Comparing DNA Sequences by word2vec Embedding.

The comparison of DNA sequences is of great significance in genomics analysis. Although the traditional multiple sequence alignment (MSA) method is popularly used for evolutionary analysis, optimally aligning k sequences becomes computationally intractable when k increases due to the intrinsic computational complexity of MSA. Despite numerous k-mer alignment-free methods being proposed, the existing k-mer alignment-free methods may not truly capture the contextual structures of the sequences. In this study, we present a novel k-mer contextual alignment-free method (called kmer2vec), in which the sequence k-mers are semantically embedded to word2vec vectors, an essential technique in natural language processing. Consequently, the method converts each DNA/RNA sequence into a point in the word2vec high-dimensional space and compares DNA sequences in the space. Because the word2vec vectors are trained from the contextual relationship of k-mers in the genomes, the method may extract valuable structural information from the sequences and reflect the relationship among them properly. The proposed method is optimized on the parameters from word2vec training and verified in the phylogenetic analysis of large whole genomes, including coronavirus and bacterial genomes. The results demonstrate the effectiveness of the method on phylogenetic tree construction and species clustering. The method running speed is much faster than that of the MSA method, especially the phylogenetic relationships constructed by the kmer2vec method are more accurate than the conventional k-mer alignment-free method. Therefore, this approach can provide new perspectives for phylogeny and evolution and make it possible to analyze large genomes. In addition, we discuss special parameterization in the k-mer word2vec embedding construction. An effective tool for rapid SARS-CoV-2 typing can also be derived when combining kmer2vec with clustering methods.